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Abstract. The possibility that the magnetic shear-flow insta- 
bility (MRI, Balbus-Hawley instability) might give rise to tur- 
bulence in a cylindric Couette flow is investigated through nu- 
merical simulations. The study is linear and the fluid flow is 
supposed to be incompressible and differentially rotating with 
a simple velocity profile with = a + b/R 2 . The simplic- 
ity of the model is counterbalanced by the fact that the study 
is fully global in all three spatial directions with boundaries 
on each side; finite diffusivities are also allowed. The investi- 
gation is also carried out for several values of the azimuthal 
wavenumber of the perturbations in order to analyse whether 
r*fS non-axisymmetric modes might be preferred, which might pro- 
fNl duce, in a non-linear extension of the study, a self-sustained 
5_i ' magnetic field. 

We find that with magnetic field the instability is easier to 
' be excited than without magnetic field. The critical Reynolds 
qq number for Pm= 1 is of order 50 independent of whether the 
t-H ' nonmagnetic flow is stable or not. Also we find that i) the mag- 
netic field strongly reduces the number of Taylor vortices, ii) 
r ~ l the angular momentum is transported outwards and iii) for fi- 
, nite cylinders a netto dynamo-alpha effect results which is neg- 
' ative for the upper part of the cylinder and which is positive for 
(f) , the lower part of the cylinder. 

■^1" For magnetic Prandtl number smaller than unity the critical 
<0 Reynolds number scales with Pm~° 65 . If this was true even 
for very small magnetic Prandtl number (e.g. liquid sodium) 
-«^. i the critical Reynolds number should reach the value 10 5 which 
i-£h ' is, however, also characteristic for nonlinear finite-amplitude 
hydrodynamic Taylor-Couette turbulence - so that we easily 



Or 



O ' have to expect the simultaneous existence of both sorts of in- 



CO 



stabilities in cold accretion disks. 
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1. Introduction 

In recent years the generation of turbulence in accretion disks 
is thought as basing on the magnetorotational instability (MRI, 
Balbus-Hawley instability) in which a differential rotation with 
outwards decreasing basic rotation rate becomes unstable un- 
der the influence of a given magnetic field (Velikhov 1959; 
Chandrasekhar 1961). After the first work (Balbus & Hawley 
1991) several other papers have studied the nonlinear evolu- 
tion of the instability (Hawley & Balbus 1991; Brandenburg 
etal. 1995; Hawley etal. 1995; Matsumoto & Tajima 1995). 
By other investigations the importance of boundaries has been 
found (Curry et al. 1994; Curry & Pudritz 1995, 1996; Kitchati- 
nov & Riidiger 1997; Kitchatinov & Mazur 1997; Riidiger et al. 

1999) . All the studies, however, were linear, a fully global non- 
linear approach being just at the beginning (e.g. Drecker etal. 

2000) . 

On the other hand, Dubrulle (1993), Richard & Zahn 
(1999) and Duschl etal. (2000) focus attention to existence of 
instabilities in hydrodynamic Taylor-Couette flow, i.e. a config- 
uration between two rotating cylinders (see also Koschmieder 
1993). There we have to distinguish two cases. For sufficiently 
slow rotation of the outer cylinder and sufficiently rapid rota- 
tion of the inner cylinder there is indeed a linear instability, 
which, however, is not interesting for astrophysical applica- 
tions such as for accretion disks. Here, however, one can ask for 
the nonlinear instabilities due to finite-amplitude disturbances. 
They indeed exist, but for Reynolds numbers exceeding those 
for linear instability by several orders of magnitudes. For a 
wide gap between both cylinders, the data discussed by Richard 
& Zahn (1999) lead to a Reynolds number for finite-amplitude 
instability of about 2 • 10 5 (see below). 

In the present paper we shall demonstrate for Taylor- 
Couette flow for plasma under the presence of an external 
field that there are always linear instabilities for rather mod- 
erate Reynolds numbers independent of whether the hydrody- 
namic Taylor-Couette flow is unstable or not. We are able to 



compute the MHD instability for magnetic Prandtl numbers of 
0.01 < Pm < 1, however, for cold accretion disks and tech- 
nical liquid metals like sodium the magnetic Prandtl numbers 
are much smaller. For a (realistic) Prandtl number of 10~ 5 an 
extrapolation is needed leading to a Reynolds number of about 
10 5 . So we have the very interesting situation that for natural 
electrically conducting fluids the critical Reynolds number for 
linear MHD-turbulence approaches the critical Reynolds num- 
ber for nonlinear hydrodynamic turbulence. If this is true then 
the theory of protoplanetary disks must be modified in an im- 
portant sense. 

We investigate here the effect of the instability in a simple 
numerical model for a cylindric Taylor-Couette flow in which a 
magnetised fluid is contained between two finite cylinders. We 
consider the permissible stationary circular flows of an incom- 
pressible fluid between two rotating coaxial cylinders (Fig. |l]). 
As we shall see, in the absence of viscosity, the class of such 
permissible flows is very wide: indeed, if denotes the angular 
velocity of rotation about the axis, then the equations of mo- 
tion allow f2 to be an arbitrary function of the distance R from 
the axis, provided the velocities in the radial and the axial di- 
rections are zero. But if viscosity is present, the class becomes 
restricted: in fact, in the absence of any transverse pressure gra- 
dient, the most general form of £1 which is allowed is 
b 



fl(R) = a 



(1) 



where a and b are two constants which are related to the angu- 
lar velocities i7j n and J7 out with which the inner and the outer 
cylinders are rotated. Thus, if R in and i? ou t (> -Rin) are the 
radii of the two cylinders, then 



b/Rl 



b/R\ 



out - 



Solving for a and b in terms of Vt- m and f2 Q ut, we have 



a = Q; 



b = n h 



#in (Wi) 



(2) 



(3) 



I — T] z 1 — l] z 

where fi = f2 out /f2 in and i) = R in / Rout- 
Flow and field perturbations are developed after the az- 
imuthal Fourier modes exp(im0). Although the flow and field 
quantities are then complex numbers only their real parts have 
an own physical meaning. 
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and 



Fig. 1. Geometry of the cylindric Couette flow 



2. Equations and model 

We solve the incompressible, dissipative, linearised MHD 
equations numerically in cylindrical coordinates for a global 
model. The basic rotation law is given in ([[]) and the rotating 
fluid may be threaded by a strictly vertical magnetic field 
B = B e z (4) 
(cf. Stone & Norman 1994). Extra radially dependent az- 
imuthal fields are also possible in this approximation (Curry & 
Pudritz 1995; Ogilvie & Pringle 1996; Terquem & Papaloizou 
1996; Papaloizou & Terquem 1997). By linearising around the 
equilibrium state and using dimensionless quantities, our equa- 
tions read 



du 



Ap 



Cn(-A-Vp) 



Ha-Tm 

Cn 



Pm Am, 



A 



Ha 2 Pm 

r 2 



Pm 



(5) 



(6) 



(7) 



— = curl(Cb (n + £)- curLB), 
at 

where p is the pressure divided by the density, u the velocity 
perturbation and B the magnetic perturbation. The normalisa- 
tions are the same as in Riidiger etal. (1999). 

In the equations A and C are the linearised advection term 
and Lorentz force given by 
q I \mu R - 2u0 
A= — I imu 4 , + 2f i u R | (8) 
m \ imuz 

and 

dB R _ dB z 
dz dr 



while 



n 



— -r~£>Z 



n 



dB^ 
dz 



f2 1 T 




(9) 



(10) 




(11) 



are the linearised electromotive force induced by the basic dif- 
ferential rotation and the flow perturbations. In the equations 
above small r's are the dimensionless radius r = R/H. Fi- 
nally, 

Ha = ^—^= (12) 

is the Hartmann number (77, p and v being the magnetic diffu- 
sivity, density and viscosity); 

Cn = — (13) 

V 

is the shear-flow parameter ('dynamo number') and Pm is the 
magnetic Prandtl number, Pm= vji). The quantity H is the 
half-thickness of the cylinder which is also the unit length for 
the non-dimensionalization. The unit time is the magnetic dif- 
fusion time t = H 2 /-q, velocities are normalized with HSl m . 
It is a time-stepping code using a global energy quenching pro- 
cedure of the form Co = Cq/(1 + E) to find the lowest 
excitation regime, where E is the total energy of the system 
(see Elstner etal. 1990). In our computations we use another 
parameter named Reynolds number, it is defined as 

r^i n _/?i n (_/? ut R\n) Ax 

KC = , (i'+j 

V 

so finally we get 

R °- PmP (15) 
with r = H/R out . 

With our definitions the nonlinear instability described by 

Richard & Zahn (1999) - using experimental data of Wendt 

(1933) and Taylor (1936) - starts for 

6 rXl-r}) 2 



Re > 5 • 10 6 



(16) 



(l + r?)3(i_ A )< 
a number which is about 2 • 10 5 for jl — and f] — 0.5 used in 
our computations below. 

A staggered mesh of 8 1 by 8 1 grid points was used to solve 
the system. In the Eqs. (§)-(0) Ha, Cq and Pm are free param- 
eters. We solved the equations in order to get a representation 
of the curves of marginal stability in the Ha-Co plane. 

We impose rigid boundary conditions for the velocity field, 
both taken on all boundary layers and set pseudo-vacuum 
boundary conditions on the magnetic field, i.e. the components 
of the magnetic flux tangential to the surface are vanishing. 
This choice also comes for simplicity in imposing the condi- 
tions for the flow, since by using a staggered mesh the grid- 
points for the velocity are not defined directly on the bound- 
aries. It also ensures the necessary vanishing of the Poynting 
flux vector in normal direction (cf. Ruzmaikin et al. 1988). It is 
generally accepted that the MRI is influenced by the presence 
of the boundaries (Curry & Pudritz 1995). It is thus of high 
importance to find out how the various boundary conditions in- 
fluence (say) the symmetry of the solutions with respect to the 
rotation axis or with respect to the equator. Local simulations 
for closed boxes can not answer such questions with their high 
relevance to observations. 
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Fig. 2. Neutral- stability lines for various azimuthal wavenumbers (m = 0, 1, 2) and rigid boundary conditions. The outer cylinder 
is at rest (left, hydrodynamic instability exists) or it rotates with 33% of the rotation rate of the inner cylinder (right, hydrodynamic 
instability does not exist). Pm= 1, t) = 0.5 



3. Results 

3.1. The neutral- stability lines 

The neutral-stability lines for to < 3 in the Re-Ha plane for 
Pm = 1 and various fi are given in Fig. ||. The curves have an 
intersection with the Re-axis; for the line to = 0, when Ha = 
0, we can get Re = 69.75. This case is well-known for the 
nonmagnetic flow, from Chandrasekhar (1961) we find Re = 
68.3. Both values are close to each other and our numerical 
simulations are consistent with the experiment. There is also 
a minimum at Re = 49.15 for the line of m = 0. We can 
conclude that with magnetic field the instability is easier to be 
excited than without magnetic field. For stronger magnetic field 
the critical Reynolds number again increases, so that as known 
the instability is suppressed by high magnetic field amplitudes. 

In Fig. |] also the case of a flat rotation law is discussed 
where no hydrodynamical instability exists (p, — 0.33). There 
is now no eigenvalue along the vertical axis. But with magnetic 
field the minimum of the Reynolds number is hardly changed 
(Re = 46.9) compared to the former case. It can be stressed 
that for Hartmann numbers of order 100 there is always a crit- 
ical Reynolds number of order 50 above which we have insta- 
bility - independent of whether or whether not a hydrodynamic 
instability exists. 

The results might be applied to two different cases. The first 
one is the solar tachocline with its values of v ~ 10 cm 2 /s and 
7/ ~ 10 3 cm 2 /s (Pm= 10~ 2 , see Riidiger & Kitchatinov 1996). 
The Reynolds numbers are much higher than 10 2 . It follows 
Bq = 10 _9 Ha/e Gauss with e = H/R for the characteristic 
magnetic field strength. With Ha ~ 10 2 and e ps 0.1 one finds 
B ~ 1 /iGauss. 

On the other hand, an experiment with fluid sodium with 
Pm= 10" 5 (y = 10~ 2 cm 2 /s, r\ = 10 3 cm 2 /s) leads to B « 
10" 1 Ha for H ~ 1 m, i.e. B ~ 10 Gauss. With R ln and AR 
of 10 cm each, it is Re~ 10 5 / with the rotation frequency /, 



FLOW 



1.01 1 1 1 _ i 1 - - -■ — 1 - - j - 




0.10 0.12 0.14 0.16 0.18 0.20 
RADIUS 



Fig. 3. The pattern for the the azimuthal flow without magnetic 
field (Ha=0) but for resting outer cylinder (ft — 0). fj = 0.5, 
Pm=l,m = 

so that it is no problem to exceed critical Reynolds numbers of 
order 100. 

3.2. The eigenfunctions 

It is known that the cell size in z for the linear nonmagnetic 
Taylor-Couette flow is practically the same as it is in r - so 
we can find 10 cells in Fig. |^ for Ha= 0. With magnetic field, 
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however, only 5 cells develop antisymmetric with respect to 
the equator (Fig. Q). Obviously, the value of the vertical wave 
number strongly decreases for increasing magnetic field (see 
Chandrasekhar 1961). 

In Fig. |] the induced toroidal magnetic field is also shown. 
It is an equatorially symmetric ('S') field. Correspondingly, the 
associated zonal kinetic flow has an antisymmetric structure 
with respect to the equator. The induced fields and flows here 
are concentrated to the inner boundary. 

3.3. Viscosity-alpha and dynamo-alpha 

Although our study is linear, it is worth analysing the efficiency 
of the instability for the angular momentum transport. The 
effective viscosity that should bring the angular momentum 
transport is generally parameterised through the quantity ass 
(Shakura & Sunyaev 1973), defined in terms of the Reynolds 
and Maxwell stress tensors by 



UrU^ 



BrB^ 



an 



with the normalisation 

v T = a ss H 2 fl, (17) 
where vt is called the eddy viscosity and ass is the viscosity- 
alpha. As we know, the ass is n °t trivially positive, there 
are also hydrodynamic simulations with negative correlations 
(u R u<j,). 

In dimensionless units 

Ha 2 Pm \ 

—^-BrBA . (18) 

In our linear approach the amplitudes of the quantities are un- 
known as they are free of a common positive or negative factor. 
Hence, the sign of quadratical terms can be computed as well 
as ratios of fluctuating quantities. 



ass oc < uru^ 




0.14 0.16 
RADIUS 



Fig. 5. The viscosity-alpha, ass, is always positive, i.e. the an- 
gular momentum transport is outwards, even for the axisym- 
metric cells (m = 0). Ha= 200, ft = 0.33, fj = 0.5 



<c -0.010 




Fig. 6. The nonlinear quantity d non-averaged as z -profile for 
m = and at a given radius. Pm= 0.1, Ha= 200, fi = 0.33, 
fj = 0.5 

Pm = 0.1 Ha = 200 m = ^ = 0.33 rj = 0.5 



B') which is due to the external field, 8 Z = a zz B z . Apart 
from positive factors we find 

a zz oc (urB^ - u^Br) (19) 
where the averaging procedure must distinguish between upper 
and lower part of the cylinder. As known the alpha-effect is an- 
tisymmetric with respect to the equator, the question is whether 
it is positive at the upper part and negative at the lower part or 
v.v. 

The results are presented in Figs. || and ^. Figure || for a 
prescribed radius displays the vertical profile of d = u^B^ — 
u<j,Bji for m — which has to be averaged to get the alpha- 
effect. The result is typical for the nature of this effect. The 
lower part of a cell has a negative (positive) d in the upper 
(lower) part of the cylinder. The upper part of the cell is op- 
posite. Due to the vertical stratification, however, the d of the 
lower part of a cell dominates the d of the upper part of a cell. 
The total alpha for each cell is negative (positive) for the up- 
per (lower) part of the cylinder. The alpha-effect only exists if 
some kind of stratification exists. Figure [7| gives the resulting 
alpha-effect after averaging d about the entire upper part of the 
cylinder or the lower part, resp. Note that our computations 
lead to a negative (positive) total alpha-effect for the upper 




Fig. 7. The dynamo-alpha is negative in the upper hemisphere 
(solid) and positive in the lower hemisphere (dashed) 



(lower) part of the cylinder. This result seems to be character- 
istic for the magneto-rotational instability (Brandenburg 1998) 
although we have been here no negative buoyancy. Obviously, 
for an infinite cylinder (or a model with periodic boundary con- 
ditions in z) there is no vertical stratification and, therefore, no 
alpha-effect. 
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Fig. 4. The eigenfunction for the toroidal magnetic field (left) and the azimuthal flow (right) in the gap between the cylinders. Ha 

= 200, Pm = 0.1, m = 



4. The magnetic Prandtl number dependence 

In order to approach real plasma we have to switch to mag- 
netic Prandtl number smaller than unity. The solar gas beneath 
the convection zone approaches Prandtl numbers of order 10 -2 
while liquid sodium due to its small viscosity with 10~ 5 is even 
more extreme. 

Simulations with Pm= 0.1 for a flat basic rotation law are 
given in detail in Fig. ||. The lowest eigenvalue is now about 
171. Some computations with Pm= 0.01 are also done. Figure 
^ yields a slope such as 

Re crit oc Pm" 65 , (20) 
which we have only justified, however, down to Pm= 0.01. If 
it would also be true for magnetic Prandtl numbers of order 
10~ 5 , the critical Reynolds number results as 9-10 4 . This value 
is very close to the Reynolds number which we know from 
( 16 ) for nonlinear but pure hydrodynamical instability. It might 
be the case that in the limit of very small magnetic Reynolds 
number, i.e. for very cold accretion disks we have the realisa- 
tion of a mixed instability: linear magneto-rotational instability 
and nonlinear Taylor-Couette flow turbulence. 

5. Discussion 

We have considered the Taylor-Couette flow experiment of a 
electrically conducting fluid under the presence of an external 
field. It is shown that for both flat or steep rotation law the 
MHD Taylor-Couette flow is always unstable. In contrast to the 
hydrodynamical regime the critical Reynolds number for insta- 



m 1000 - 



500 - 




\m = 2 



400 600 
HARTMANN NUMBER 



Fig. 8. Neutral-stability lines for various azimuthal wavenum- 
bers (to = 0,1,2) and rigid boundary conditions with Pm 
= 0.1. Hydrodynamic instability does not exist. The minimum 
Reynolds number is about 176. ft = 0.33, i) = 0.5 

bility does hardly depend on the rotation of the outer cylinder. 
The dependence of the critical Reynolds number on the mag- 
netic Prandtl number, however, is much stronger, see Eq. (|20|). 
It is probably not too strong to avoid experimental work with 
liquid sodium. 

For Pm= 10~ 5 (liquid sodium) the estimated Reynolds 
number is of order 10 5 , which is also known from experi- 
ments concerning the nonlinear instability of Taylor-Couette 
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-2 -1 

og Pm 



Fig. 9. Critical Reynolds numbers for various magnetic Prandtl 
numbers. The Reynolds number increases for decreasing 
Prandtl number, m = 

flow. So we have the puzzling situation that the Reynolds num- 
ber for two very different types of instability could be very sim- 
ilar, i.e. for nonlinear hydrodynamic instability and linear mag- 
netohydrodynamic instability. If true, such a result will have 
strong implications for the simulations of protoplanetary disks, 
which due to the low temperatures also have very low magnetic 
Prandtl numbers. 

Of particular interest for experiments also should be the ex- 
istence of a turbulence electromotive force. Although its mag- 
nitude here remains open (as the theory is linear) we can find 
the sign of the bilinear expression. The EMF divided by the 
uniform external field the dynamo-alpha is found for the con- 
tainer. It is the zz-component of the alpha-tensor as only a ver- 
tical magnetic field is applied. The a zz proves to be negative 
(positive) in the upper (lower) part of the cylinder. It only ex- 
ists due to the finiteness of the container in vertical direction. 
One can understand this finding with the following. In order to 
construct the pseudo-tensor a one needs a scalar g ■ ft with a 
characteristic direction g. In an vertically infinite cylinder the 
only characteristic direction g is the radial one, but g ■ fi = 
in this case. Only for cylinders which are bounded in the verti- 
cal g || f2 acts as an characteristic direction leading to a finite 
alpha-effect with an antisymmetry with respect to the equator. 
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